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In the framework of the Ginzburg-Landau approach, we present a self-consistent theory of specific 
soliton states in mesoscopic (thin-walled) two-band-superconducting cylinders in external parallel 
magnetic fields. Such states arise in the presence of " Josephson-type" interband coupling, when 
phase winding numbers are different for each component of the superconducting order parameter. 
£Sj ■ We evaluate the Gibbs free energy of the sysyem up to second-order terms in a certain dimensionless 

_ ■ parameter e rn <C 1, where C m and Ck are the magnetic and kinetic inductance, respectively. 

We derive the complete set of exact soliton solutions. These solutions are thoroughly analyzed 
from the viewpoint of both local and global (thermodynamic) stability. In particular, we show 
that rotational-symmetry-breaking caused by the formation of solitons gives rise to a zero- frequency 
rotational mode. Although soliton states prove to be thermodynamically metastable, the minimal 
energy gap between the lowest-lying single-soliton states and thermodynamically stable zero-soliton 
states can be much smaller than the magnetic Gibbs free energy of the latter states, provided that 
intraband "penetration depths" differ substantially and interband coupling is weak. The results of 
our investigation may apply to a wide class of mesoscopic doubly-connected structures exhibiting 
. two-band superconductivity. 
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C3 , I. INTRODUCTION 



The subject of this paper is a self-consistent theory of specific soliton states that were originally predicted in Refi 
and reportedly observed experimentally. 2 Without any doubt, these states can be regarded as a hallmark of two-band 
q , superconductivity in mesoscopic doubly-connected samples. 

Indeed, owing to the emergence of additional degrees of freedom of the order parameter, the nomenclature of topo- 
L " "™ ', logical objects in multiband superconductors is much richer than that in conventional single-band superconductors. 
I ■ In particular, Ginzburg-Landau equations describing two-band superconductivity in bulk samples admit topologically 
stable solutions (with one-dimensional singularities of the order parameter) that can be interpreted as vortices carrying 
■^j- ■ fractional magnetic flux^ In the absence of any interband coupling, these vortices are accompanied by a circulating 
neutral superfow associated with gradients of the interband phase difference. In the presence of "Josephson-type" 
interband coupling, the neutral superflow generates static solitons of the sine-Gordon type. In contrast to tradi- 
tional Abrikosov vortices in type-II superconductors, the energy per unit length of these composite topological defects 
diverges at spatial infinity: hence they are thermodynamically metastable and difficult to create in bulk samples. 

However, solitons of the interband phase difference can exist by themselves in doubly-connected mesoscopic samples, 
when the formation of any magnetic vortices in the volume of the superconductor is prohibited energetically. 1 Moreover, 
soliton states in this case can be induced by an externally applied magnetic field, which makes them a convenient 
object of investigation. Thus, experimental studies^ of the magnetic response of mesoscopic two-band superconducting 
rings reveal certain nontrivial features that, according to the authors of Ref^, can be attributed to the creation of 
metastable soliton states. 

Our research is largely motivated by the absence in current literature of any quantitative theoretical analysis of 
this pronounced feature of two-band superconductivity. (Unfortunately, the arguments of Ref.i and of the recent 
publications^ are mostly heuristic by nature.) Mathematically, the approach of this paper is based on a Ginzburg- 
Landau-type theory, which is a commonplace in theoretical studies of topological defects in two-band superconductors: 
see the next section. This means, of course, that we are restricted to the temperature range 

T c -T 

where T c is the critical temperature of the superconducting transition. 

As to the physical object, we consider a two-band superconductor in the form of a straight, circular thin-walled 
cylinder, whose symmetry axis is the z axis of cylindrical coordinates (r, (p, z) (see Fig. [T]). The constant external 
magnetic field H is applied along the symmetry axis: H = (0, 0, H > 0). The length of the generatrix of the wall of 
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FIG. 1: The geometry of the problem (schematically). The parameters L, R and d obey conditions ([2])-([5]). 
the cylinder satisfies the condition 

L»R=^, (1) 

which allows us to neglect end effects. The wall thickness d = — r\ and the average radius of the cylinder R satisfy 
the following conditions: 

d«min{a,6}, (2) 
fl> max {£!,&}, (3) 
R > A, (4) 

where £i and £2 are the "coherence lengths" in bands 1 and 2, respectively, and A is the weak- field penetration depth. 
Exact definitions of £1, £2, and A will be given in the next section; however, the role of conditions ©- © should be 
explained right now. Thus, condition ([2]) precludes the formation of any magnetic vortices in the wall of the cylinder. 
Condition (UJ) mainly simplifies mathematical consideration. In contrast, a combination of conditions Q and (J5j is 
of primary importance: taken together, these two conditions guarantee that self-induced magnetic fields are small 
and can be treated perturbatively. (This fact justifies the definition "mesoscopic cylinders" in the title of the paper.) 
Moreover, as will be shown, the dimensionless quantity e serves as a natural expansion parameter of the Gibbs free 
energy. In order to carry out a rigorous analysis of thermodynamic stability of soliton states, we will have to evaluate 
the Gibbs free energy exactly up to small terms of order e 2 , which implies the necessity of self-consistent evaluation 
of the vector potential up to first-order terms in e. 

We conclude the formulation of the problem by specifying boundary conditions for soliton states. Consider a 
two-component superconducting order parameter ^ — (^1,^2)) where ipi = \tpi \ e 1 ^ 1 and ip2 — IV^I e I( ^ 2 . The double- 
connectedness of the cylinder is accounted for by the condition of single- valuedness of In particular, this condition 
requires that 

IV>l|L=o = l^l|L=27T' l^2|L =0 = I^IL= 2w - ( 6 ) 
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As to the phases <f>x and (f> 2 , the requirement is as follows: 

j> V0i • dl = 2-imx, <j)V(f) 2 ■ dl = 2nn 2 , m, 2 = Q,±l,±2, . . . , (7) 

where T is an arbitrary closed continuous contour that lies inside the wall of the cylinder and encircles the opening. 
It should be emphasized that there are no a priori reasons for setting ri\ = n 2 & As in the case of fractional magnetic 
vortices in bulk two-band superconductors^ nontrivial topological states arise when n\ ^ n 2 . In the presence of 
interband coupling, they are of the soliton type. 

In Section II, we introduce the Gibbs free-energy functional of the system and analyze its basic properties. In 
Section III, we derive a self-consistent expression for the electromagnetic Gibbs free energy. Soliton solutions are 
derived and thoroughly discussed in Section IV. Finally, in Section V, we summarize the obtained results and make 
several concluding remarks. Appendices A and B contain details of some mathematical calculations skipped over in 
the main text. In Appendix C, we present several particular examples of soliton solutions. 



II. GIBBS FREE-ENERGY FUNCTIONAL 



We begin by defining the Gibbs free-energy functional of the system. In view of complete homogeneity along the z 
axis and with the normal-state Gibbs free energy being subtracted, it takes the following form: 

G [*,**, A; H] =L I d 2 v a x |^| 2 + a 2 |0 2 | 2 + ^ |Vi| 4 + ^ |0 2 | 4 
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Here, the coefficients (3x and j3 2 are positive constants, whereas ot\ and a 2 are temperature-dependent: 

a\ = ai (T) = a x (T - Ti) , ax = a 2 (T) = a 2 (T - T 2 ) , ax, a 2l Tx, T 2 > 0. 
Moreover, the latter coefficients enter the definitions of the coherence lengths £i and £2: 

h h 
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V2miai ' "*° s/2m 2 a 2 
The electron charge in © is e < 0; the total vector potential A defines the local magnetic field h: 

h = VxA, h=(0,0,/i), h = h(r). 



(8) 



(9) 



(10) 



(11) 



The parameter of interband coupling, 7, may have either sign. Two-dimensional integration in the plane (r, ip) is 
carried out over the cross-section of the superconductor (Eg) in the square-bracketed terms, and over the cross- 
sections of the superconductor and of the opening (Ss + So) in the last (magnetic) term. 

A microscopic derivation of free-energy functionals of the type © was given in Refji for the case of clean two- 
band superconductors in the limit of small interband coupling. Free-energy functionals of this type are employed in 
theoretical studies of different aspects of two-band superconductivity, such as, e.g., topological defects current- 
carrying states^ the Little-Parks effect^ surface energy^ etc. It should be additionally noted that, for our specific 
geometry and 7 > 0, the functional © also applies to the description of a composite system consisting of two thin 
coaxial cylindrical films of single-band superconductors, Josephson coupled via a parallel insulating layer, 12 which is 
exactly the experimental set-up of ReL 2 -. 

To obtain the actual (observable) Gibbs free energy, one has to minimize © with respect to "J, ^* and A under 
appropriate boundary conditions; however, substantial simplifications can be made already at this stage. First, we 
notice that, by the symmetry of the problem, the amplitudes |-0i| and |^2| do not depend on ip. Moreover, they 
cannot depend on r, either. Indeed, by virtue of condition ([2]), any radial variations of |-0i| and |-02| would give 

rise to free-energy terms that are by the factors % S> 1 and % 3> 1 larger than the first and the second terms in 
([5]), respectively, which is energetically unfavorable. 13 As a result, in equilibrium, the magnitudes \ipx\ and |^2| are 
functions of T and H only. 
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Consider now the kinetic-energy terms [the first two terms in the second line of JSJ]. The ratio of these terms to 

the first and the second terms in (jH]), respectively, is at most of order (~*o"J an< ^ W (~*o~) ' wnere ®h is the 
external flux, and 

— W 

is the flux quantum. Owing to condition ([3j> , for sufficiently weak external fields, ^$ (^f) ' W (^o~) ^' (C om - 
pare with the consideration of flux quantization in singe-band-superconducting cylinders. 14 ) In this field range, we 
can set |^i| = IV'ilo an d \ip2\ = IV^Io: where IV'ilo an d |"02 lo satisfy the equilibrium conditions for an unperturbed 
two-band superconductor: 

ai|M>+ftMo-|7llM> = 0, 

^I^Io + ^I^Io-ItII^iIo-o. (13) 

One can readily obtain a good approximate solution to (|13p . 15 However, it is of no interest in the context of the soliton 
problem. We only note that the critical temperature, derived from ([T3|) . is T c = | Ti + T 2 + >/ (Ti — T2) 



In light of these simplifications, it is reasonable to consider the weak-field penetration depth 10 
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A=t^oJ ■ .2 .2 ( 14 ) 
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and to define intraband "penetration depths 

Ai ~v^mW 1 + 2 " ■ (15) 

(For the above-mentioned composite, Josephson-coupled system, the quantities Ai and A 2 have direct physical mean- 
ing.) We also introduce new, functionally independent phase variables fi> and ^^ilii 

</> = (pi-<P2, (16) 

X = cicpi + c 2 (p2] ci = (AAf 1 ) 2 , ^^(AA^ 1 ) 2 , ci+c 2 = l. (17) 

Using definitions (TT^j) and (fT4]) - ([TT)) . we obtain the reduced Gibbs free-energy functional in the following form: 

G [ct>, X , A; H] = F S0 + G em [ X , A; H] + F sol [4>] . (18) 

Here, the first term is the free energy of the unperturbed superconducting cylinder) 1 ^ 

f so = v s L \Hl + «2 \Hl + y l^ilo + y l^lo - 2 M IV'ilo l^lo) ; ( 19 ) 

V s = 27ri?L<i. 

The second term is the electromagnetic Gibbs free-energy functional: 

0_ k , A;H] = J^/r (v x - |a)' + ± f^f, (h - H)> , (20) 

with the first term on the right-hand side of (|20[) being the kinetic-energy functional of the supercurrent. Finally, the 
last term in (fl8l) is 



S2 , , 

(21) 



2 

(V0) + ^ (1 — sgn 7 cos i 



/2 _ *o 1 c i c 2 



32tt3 A 2 MIViloWo 

where sgn x is the sign function. The term (f2"Tj) should be interpreted as the soliton self-energy functional. Indeed, 
when n\ — n% in ([7]), we have^^ either (/> = mod 2ir (for 7 > 0) or = 7r mod 27r (for 7 < 0), and this term vanishes 
identically. 

Our task now is to minimize (|18l) with respect to 0, x an d A. As the phase variable <j> is not coupled to the vector 
potential A, this procedure can be performed in two separate steps. 
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III. ELECTROMAGNETIC GIBBS FREE ENERGY 



The minimization of the electromagnetic functional (|20[) reduces to evaluation of the stationarity condition SG em = 
0, or 

^.=0, ^L =0 . (22) 
oA ox 

Indeed, in view of quadratic nature of (|20p . solutions to (|22p are automatically minimizers of this functional (i.e., 
5 2 G em > at these solutions). 

Variation with respect to A yields Ampere's law 

VxVxA = 0, re (0,71); (23) 

4-7T 

VxVxA = — j, re (run), (24) 
c 



j = -7^(^V X + A) ,25) 



with 

~ 4ttA 2 V 2?r 

being the supercurrent density [j = (0,j, 0) by symmetry], and the boundary condition 

h| r=r2 = V x A| r=ra = H. (26) 

(This boundary condition should, of course, be complemented by the conditions of continuity of A and h at r = r 2 and 
the condition of regularity of A at the origin.) Variation with respect to \, under the condition of single- valuedness 
of variations Sx, just yields the current-conservation law 

Vj = (27) 

and the single- valuedness condition 

jVo^W- (28) 

[This boundary condition should be complemented by a condition on x resulting from 0.] 

The problem of finding A and x is still sub-definite, because we have not so far fixed the gauge. As the z component 
of the vector potential drops out of the definition of h [see (|lip]. it is equal to an arbitrary constant, and we set A z = 0. 
The r component of the vector potential can be eliminated by the gauge transformation 

A^A-V/ A r (r',<p)dr', X -> X + \ A r (r' , <p) dr' . 
Jo ^oJo 

In this particular gauge, 

A=(0,A,0), A = A(r); (29) 

h=~(rA), (30) 
r ar 

and x does not depend on r (j r = 0). Using ([7]), (fT7|) . (I2T1) and (j2"5)l . we arrive at a well-posed boundary- value 
problem, 



fx 



0, ^€(0,2^): 



X (2tt) = x (0) + 2vr (ma + n 2 c 2 ) , ^ (2tt) = ^ (0) , 



whose solution is 



X (</?) = ("lCl + ^2C 2 ) (y 5 + <A), 



(31) 
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with ifQ being an arbitrary constant. 

The boundary-value problem for the vector potential now takes the form 



d 
dr 



r dr 



d 
dr 

(rA) 



1 d , ^ 
-— (rA) 
r dr 



0, re(0,n); 



1 

A2 



A J. $0 < \ 



re (n,r 2 ) 



q(n 1 ,n 2 ) = n Y c x + n 2 c 2l 
1 d 



|A|| P=0 < oo, A\ r=ri _ = A\ r=ri+0 , ~ (rA) 



r— n — 



1 d 

r dr 



(rA) 



r=r i+0 



-— M 

r dr 



= H. 



(32) 



This boundary-value problem admits an exact solution: it is presented in Appendix A. However, to obtain a second- 
order expansion of the electromagnetic Gibbs free energy in terms of the small parameter e (see Introduction), we 
need only first-order expansions of A and h. They are as follows: 



v r 
A(r) = -H-- 
W 2 2 

_r n 



q (711,712) + 



$0 f \ _l ft' 
— 2? (7ll,«2) + -n 
7rr^ 



e, re [0,n] ; 
e, re (ri,r 2 ] ; 



(33) 



ft(r) 



if 



$0 

— o<7 (ni,n 2 ) 
nrf 

r 2 -r f $ 



if 



£, re [0,ri] 



r 2 - ri 



q(n 1 ,n 2 ) + H 



r e (ri,r 2 ] 



(34) 



[The fact that expressions (|3"3"]l and on the interval (ri,r 2 ) are not related to each other by equation (|3T)|) should 
not cause any confusion: to ensure the fulfillment of (|30p . we would have to continue the expansion of (|A4[) up to 
small terms of order and e-j|.] 

Bearing in mind that in integral physical quantities any difference between r l7 r 2 and i? should be neglected [sec 

(|A8j) ]. by use of equations ([34]) and ([24| we immediately obtain expressions for the total flux $ = / hd 2 r and the 



(35) 



total supercurrent J = L / jdr: 

$ = + $ H = nR 2 H, $i = -[$ <7(rti,n 2 ) + $ J F/]e; 



c 4^ 2 i? 2 , . 

J =7-**, £ m = — p- ■ (36) 

Here, and are the external and self-induced flux, respectively; C m is the magnetic inductance of the cylinder 
(or self-inductance) 

On substitution of relations (f3~Tj). (|3"3"|) and (|34p into (j2T))) . we get the electromagnetic Gibbs free energy of the state 
parameterized by topological numbers rii and n 2 : 

r , rrs [$og(ni,n 2 ) + $ g ] 2 [$og(rti,n 2 ) + ^ g ] 2 a 
G em (ni,n 2 ; H) = — e (1 - 2s) + — s 

= ^[q(n 1 ,n 2 ) + f H ] 2 (l-s), f H = ^-. (37) 

(Note that the first term in the first line of this equation is the kinetic energy of the supercurrent, whereas the second 
term in the same line is the magnetic Gibbs free energy.) 
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To facilitate an analysis of soliton states, we should transform (|37|) to a more convenient form. First, instead of 
parameterization by n\ and n 2 , we introduce parameterization by n\ and the soliton number n = n\ — n 2 . From now 
on, we assume (without any loss of generality) that < c% < c\ < 1. Relation (|37| is rewritten as follows: 

Gem (ni,n; f H ) = {f H + m - nc 2 f (1 - e) . (38) 

From a thermodynamic point of view, of interest is the minimum of (|38[) for given |n| = 0, 1, 2, . . ..Therefore, relation 
(|38p should be minimized with respect to n\ and sng n. With this in mind, we introduce two discontinuous functions, 
a step function m (x) and a periodic function 9 (x) , via the definitions 

( s / [as], 0< {x} < i; , . 

m i = , ' i F t 1 (39) 
w \ [ar] + 1, | < {4 < 1, v ; 

and 

where [x] and {x} are the integer and fractional parts of x, respectively. Given that fjj and nc 2 can now be represented 
as Jh = m (/h) + $ (/ff) an d =sgn n [m (|n| C2) + (|n| c 2 )\, respectively, the result of the minimization is 

Gem(\n\;f H )= min G em (m, n; /#) 

ni,sgn n 

= G CT „ (-m (/jj) + m (|n| c 2 ) sgn 6 (f H ) sgn (|ra| c 2 ) , \n\ sgn 6* (/ ff ) sgn (\n\ c 2 ) ; f H ) 
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{\8(f H )\-\9(\n\c 2 )\] 2 (l-e). (41) 



IV. SOLITON STATES 
A. Soliton self-energy 

The variation of (|2ip with respect to cj>, under the condition of single-valuedness of variations 
two-dimensional sine- Gordon equation in polar coordinates, 



1 d 2 4> Id 



r 2 dp 2 

and the boundary conditions 



r dr V dr 



r dip 



sgn 7 

— 2— sin^, re(n,r 2 ), <^e(0,27r). 



= 0. 



1 d<£ 




dcj) 


d(j> 




) 


dr 





yields a static 
(42) 

(43) 



[These boundary conditions should, of course, be complemented by a condition on <f> resulting from ([7)1 .] 
However, equation (|42[) . in its exact form, by far exceeds the accuracy of our calculations in the previous section [see 

expressions (|33|) . (j34|) and relations (|A8|) ]. Discarding in (|42|) terms of order ^ and we arrive at a two-dimensional 
sine-Gordon equation in "Cartesian coordinates": 



d 2 <j) ^ d 2 4> i? 2 sgn 7 . 
dip 2 dp 2 I 2 



P =R- 



(44) 



Solutions to (|44[) . minimizing the functional (|21[) . should not depend on p for symmetry reasons [which, of course, is 
compatible with boundary conditions (I43[) ]. 

Thus, the phase 4> = 4> (</?) satisfies the following boundary-value problem: 



d 2 4> R 2 sgn 7 . , „ . 

= ^-^sm(/>, pe{0,2n); 

^(2^) = 0(0) + 2,™ (n = ±l,±2,...), ^(2tt) = ^(0) 



(45) 
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The solution of (|45l) is straightforward: 17 



(1 + sgn 7 )7T / nK {k n ) 
4>n (f) = o ^ 2am vP ~ Vno) , k n , (46) 



where am it is the elliptic amplitude, 18 K (k) is the complete elliptic integral of the first kind^ p n o are arbitrary 
constants, and k n (n = ±1, ±2, . . .) satisfy the equations 

nR 

\n\k n K(k n ) = —, u = ±l,±2,.... (47) 

Particular examples of solutions (j46|) that possess asymptotics in terms of elementary functions are relegated to 
Appendix C. Nonetheless, the very special class of exact elementary solutions is worth being presented here: namely, 
the non-soliton nontrivial topological solutions corresponding to the physical case of the absence of interband coupling 
(ItI = 0)- These solutions can be obtained from (|46|) by the limit procedure 

fcn->0, n = ±1,±2,..., (48) 

and they have the general form 

4> n (ip) = nip + tp . (49) 

They are necessarily minimizers of ([2T]) [i.e., at these solutions S 2 F so i > 0, because the functional (|21[) is quadratic in 
the case I7I = 0], and their self-energy is 

(n) = ]im F sol (n) = * \n\ 2 Cl c 2 . (50) 

If I7I 7^ 0, the functional (f2Tj) is non-quadratic, and we should analyze the second variation of (J2T|) in more detail. 
To this endji 7 -^ we turn to the Sturm-Liouville problem 



(ftp 

-j-2 + cos ft> n ip = nip, ip e (0, 2tt); (51) 



where fi> n is a given solution from the set (|46)). As shown in Refsi 17 ' 19 , 

i-2-k 

S2F s°l\<P =c t, n - ^0 / l^nl 2 dp, 



where u is lowest eigenvalue of the problem ([51]) . In our case, both /io and the corresponding eigenfunction ipo can 
be readily found: 

^0 = 0, Vo = const dn \ K kn (p - p n0 ) , k n 



V 77 

where dnu = d ™ - .— This means that S 2 F so i\^_^ > 0, and soliton states turn out to be indifferently stable states. 
Indeed, the zero value of no should be attributed to the existence of a zero-frequency "rotational mode" (by analogy 
with the well-known^ translational mode in quantum field theories) that restores rotational symmetry broken by the 
formation of solitons. To prove this, consider a small variation of 4> n induced by a small variation of the constant of 
integration p n o: 



1 \ if n \ r \ if n K i^n) , \ 
(p n {(f) ->• 4> n \ ip + —a J w <p n (p) + adn I (<p - p n0 ) , k 



\a\ < 1. 



From the above, we see that 6(j> n oc ip . 

Now that local stability of soliton solutions is established, we proceed with a discussion of soliton self-energy. It is 
obtained by the substitution of solutions (|4"6")l into (p?Tj) and has the form 

Fsoi (n) = ^^-c lC2 K(k n ) [2E(k n ) - (l - h») K (k n )} , (52) 
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where E (k) is the complete elliptic integral of the second kind^ 

First, we note that the constants of integration (p n0 that figure in (|46l) drop out of the right-hand side of ([52)1 , as 
they should. The self-energy does not depend on the sign of 7 and of n, either. By considering (formally) |n| as a 
continuous variable, we get 



8F sol (n) _ $g e 4 | n | R 



d\n\ 



Cr. 



7T/ 



CiC 2 E (k n ) > 0, 



which means that F so i (n) increases monotonically with an increase in \n\, as could be expected. However, in contrast 
to the case I7I = [see (|50j)]. the growth of F so i (n) is slower than |n| , because 



d 



F S oi (n) 



$ge AR 



E(k n )-(l-kj)K(k n ) ^ 
cic 2 < 0. 



Gi 



that 



C m 7T \n\ I 

dF sol (n) _ $2 e 4 |„| 2 E (k n ) [E (k n ) - (1 - k 2 n ) K (k n )] 



dk n 



-C1C2- 



fc„(l 



>o, fc n e(o,i) : 



the self-energy increases monotonically with an increase in fc„ on the whole interval (0, 1). The minimal value of (|52 
is achieved at k n — and is given by (|50[) . In view of the relation 



dk n 

Hi) 



7TR 1 



\n\ 2 l 2 E(k n ) 



< 0, 



the self-energy decreases monotonically with an increase in <G (0, 00) (for a given |n|). [In other words, F so i is an 
increasing function of the interband coupling parameter |7|: see the definition of I in (|21l) .] 



B. Thermodynamic metastability 

According to (|T5|). (|4ip and (|5^)) the minimal Gibbs free energy of soliton states with a given |n| in the field H can 
be represented as follows: 



G(\n\;f H )=F 



4 In 



Hi 



'(fH)\-\e(\n\c 2 )\] 2 (l-s) 



+ ^-C lC2 K (k n ) [2E (k n ) - (1 - *£) X (fc„)] 



(53) 



To analyze thermodynamic stability of soliton solutions, we should compare expression (|53[) for \n\ > 1 with the Gibbs 
free energy of the states with \n\ = 0: 



G(O:f H )=F o + ^\0(f H )\ 2 (l-e) 



(54) 



With this in mind, we first note that, for |n| > 1, the energy G (|n| ; fjj) increases monotonically with an increase 
in |n|: see Appendix B for a proof. [Contrary to what may seem, this fact is by no means obvious, because the 
electromagnetic term in (|53[) may decrease with an increase in \n\.] Furthermore, since expression (|50[) provides the 
greatest lower bound for soliton self-energies, we can restrict ourselves to the case \n\ — 1 and k n = 0. Bearing in 
mind that c 2 & (0, i] by assumption (see the end of Section III), we arrive at the following important inequalities: 

AG(l-J H ) = G(l;f H )-G(0;f H ) 



> AG<°» (1; f„) = lim G (1; f H ) - G (0; f H ) 



c 2 [l-c 2 e-2\0(f H )\(l-e)} > 0. 



(55) 
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FIG. 2: Gibbs free energy of several different topological states (m, n) (|n| = 0, 1, 2) for typical values of the 
parameters c\ , C2 , and y . Note double degeneracy of soliton states in the case c\ = c% = | and see the text for 

further explanations. 



The above inequalities clearly demonstrate thermodynamic metastability of soliton states and bring to light certain 
subtle physical points. In particular, 



maxAG(l;/ ff ) = AG(1; f H )\ fa=p > w^-c 2 (1 - c 2 e) , p = 0,l,2,. 



In contrast, 



2 ,.2 



minAG(l;/ ff )= AG(l;/ ff )| /H=p+ i > |^c lC2 , p = 0,1,2,..., 



(56) 



(57) 



which shows that min/ H AG(l;/#) can be much smaller than the magnetic Gibbs free energy of the zero-soliton 
states [see (j37|) ]. provided that < c 2 <C c\ < 1 [i.e., when Ai <C A2 < 00, see ([17))] and 3> 1 (weak interband 
coupling) . 

In Fig. [21 we plot the Gibbs free energy of several different topological states (ni,n) (|n| = 0, 1,2). Thermody- 
namically stable zero-soliton states are denoted by thick solid lines. Minima of the Gibbs free energy of soliton states 
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represent the soliton self-energy and occur when the self-induced flux $i compensates for the external flux i.e., 
when 

fn +ni- nc 2 = 0, 

[see relations ([55)1 and (|41[) ]. In the very special case, when c 2 = c\ = \ and n = 2nx, no flux is induced = 0), 
and minima of the soliton Gibbs free energy occur at H = 0. 



V. SUMMARY AND CONCLUSIONS 



Summarizing, we have presented (in the framework of the Ginzburg-Landau approach) a self-consistent theory 
of specific soliton states that constitute a distinctive feature of two-band superconductivity in mesoscopic multiply- 
connected samples. Although our mathematical consideration concerns the concrete geometry of Fig. 1, the final 
results can be expressed in terms of the magnetic and kinetic inductance (see Ref J^) and, therefore, should apply to 
a much wider class of structures. This allows us to make several generalizing remarks. 

As the predicted fractional magnetic vortices in bulk two-band superconductors^ the soliton states considered here 
prove to be thermodynamically metastable. However, the minimal energy gap between the lowest-lying single-soliton 
states and thermodynamically stable zero-soliton states can be much smaller than the magnetic Gibbs free energy 
of the latter states, provided that the intraband "penetration depths" (fT5|) differ substantially and the interband 
coupling is weak. (In order to establish this important physical fact, we had to evaluate self-consistantly the vector 
potential. The results of this evaluation may be of interest in themselves.) 

Our consideration encompasses in a natural way the case of superconducting Josephson-coupled bilayer structures 
studied experimentally in Ref. 2 . Our conclusion that the self-energy of soliton states increases monotonically with an 
increase in the strength of interband coupling qualitatively agrees with the observations reported therein. 

Furthermore, as a particular limit, our consideration contains the case of zero interband coupling. In view of the 
recently discussed possibility of independent superconductivity of electrons and protons in a liquid metallic state of 
hydrogen^ some of our results may find application in this situation as well. 

Finally, the exact soliton solutions derived in this paper should be compared with the exact soliton solutions 
representing equilibrium Josephson vortices in a superconducting tunnel junction ! 17 ' 19 In particular, Josephson vortices 
are pinned in their equilibrium positions owing to interaction with the edges of the junction. In contrast, in rotationally 
symmetric doubly-connected two-band superconductors, soliton positions are not fixed, which gives rise to a specific 
zero-frequency rotational mode. However, any defects that break rotational symmetry must cause soliton pinning. 
The effect of this pinning requires a separate discussion. In conclusion, we hope that our paper will stimulate further 
experimental and theoretical studies of the intriguing phenomenon of soliton states in two-band superconductors. 

Appendix A: Exact solution of the boundary-value problem for the vector potential 

The exact solution to (f3"2"j) has the following form: 

A(r) = -^, r€[0,n]; 
= -^(n 1) n 2 ) + C 1 / 1 (^)+(7 2 ^ 1 (^), re(ri,r 2 ]; (Al) 



h (r) = ho, r £ [0,ri] ; 



(A2) 



2\Hf 2 ft, n- 



<M(ni,n 2 )Aftft) 



nh ft, + 2A/ 2 ft, rx) nn [n/i ft, ^) + 2A/ 2 ft, £)] ' 



(A3) 



XH [nKp ft) + 2XKi ft)] A$ og (ni,n 2 )iroft) 
nh ft, + 2A/ 2 ft, m [n/x ft, a) + 2A/ 2 ft, %)] 
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r _XH [ ri I ft) - 2Xh ft)] , A$o9(ni,n 2 )7 (f) 



nh ft, r ±) + 2A/ 2 ft, 51) Trn [n/x ft, ^) + 2A/ 2 ft, £)] ' 

fi (x, V) = h (x) K (y) - I (y) K (x) , f 2 (x, y) = I (x) K x (y) + I x (y) K (x) . 

Here, l v (x) and K v (x) are modified Bcsscl functions of order v = 0, 1^ ho is the constant magnetic field in the 
opening. Expressions (|Al|) - (|A3|) are greatly simplified under condition (j4]): 

A(r) = r -^, rG[0,n]; 

$ o9 K, n 2 ) sinh x + fr ( cosh x - V^r cosh ^ 



2nr sinh j + ^ cosh j 



n-cosh 1^1 + 2A S i n h E=a 

-A£M- — £ 7^, re(n,rj]; (A4) 

V r sinhf + ^coshf V J V ; 

h(r) = ho, re[0,n]; 



$0Q , (ni,n 2 ) sinh 



'r 2 — r 
A 



Trriy 7 ??! sinh j + 7^ cosh ^ 



H-Sinh I-fl + 2A CQSn r^rj. 

+ H \- r^r-Tx re(n,r a ]; (A5) 

V r sinh j + 77 cosh | 

_ *og(ni,«2) i h d + 2A /7I g 

° sinh 4 + — cosh 4 ' 1 ' 

Expressions (|A4|) - (|A6|) should be compared with analogous expressions for a single-band-superconducting cylinder.— 
As can be easily seen, the cylinder exhibits a considerable Meissner effect under the conditions 

d -, dR , . 

- x «l, (A7) 

In contrast, in the opposite case, when condition is fulfilled, the Meissner effect is small, and expressions (|A4[1 - 
(|A6[) can be readily expanded up to first-order terms in e. Taking into account a hierarchy of the small parameters 
of the problem, 

d d d _ X 2 d _ X 

7 1 ~T 2 ~li = 2£ lp , x = 2s r ; 

d d X , . , 

R < A <<: £ ' R <<C ' ( ^ 

we arrive at the first-order expressions (|33l) and (p4|) . 

Appendix B: A proof of the inequality G (|n| ; /h )|i„i >1 > G(1;/h) 
Consider the expression 

— 2jAG (M ; /h) = ^7 [G (|n| ; / fl ) - G (0; / H )] 

= -\0 (\n\ 02)| [2 |0 (/ H )| - |0 (|n| c 2 )|] (1 - e) + ^c lC2 ^ (fc„) [2B (fc„) - (l - k 2 n ) K (k n )} (Bl) 

that follows directly from Our task is to prove that the right-hand side of (|B1[) for |n| > 1 is larger than for 

|n| = 1. Given that c 2 € (0, |] by assumption (see the end of Section III), it is sufficient to provide a proof for 
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1 < n < 



1 

2c 2 



1, where 



l 

2c 2 



is the integer part of Indeed, the first term on the right-hand side of (|B1[) 



satisfies the inequality 



whereas 



\\9(\n\c 2 )\[2\e(f H )\-\9(\n\c 2 )\](l-e)\<~(l- S ), 



AK(k v 



[2E (k n ) - {I - k 2 n ) K (k n )} > 1. 



Therefore, for |n| > 



l 

2c 2 



1, when |n| C2 > s, any possible decrease in the first term on the right-hand side of (|B1 
due to an increase in |n| cannot compensate for an incurred increase in the second term. 
For 1 < Inl < 



l 

2c 2 



1, there holds the relation \n\c2 < h, and (|B1I) becomes 
2C r 



AG(|n|;/ H )| w>x 



|n| c 2 [|n| (1 - c 2 e) - 2 |0 (1 - e)] + cic 2 
, 2£ 



4if (fc„ 



[2£ (£;„)- (1 - A" (fc n )] -1 



= n 



AG(l;/ ff )+ C2 |n|(|n|-l) 



(l-c 2 e) 
>^AG(1;/k) 



f£l^l[2£(fc„)-(l-^)if(M] 



Cl 



which was to be proved. 



Appendix C: Particular examples of soliton solutions 



In Fig. 02 we present several different soliton solutions obtained numerically. (For greater clarity, we plot the 
derivatives ^ L .) However, in two limiting cases soliton solutions possess asymptotics in terms of elementary functions. 
Thus, for f < 1, we haveJI^ 



, . (I + sgn 7) 7r , , R 2 . 

VP) ~ 2 hn((p-(p ) + sin [n {ip - tp )\ 



(Cl) 



The self-energy of soliton solutions (|C1|) is 



n\ cic 2 1 



2R 2 
\n\ 2 P 



(C2) 



[Notice that expression (|CJ2[) clearly illustrates the general features of the self-energy of soliton solutions established 
in Section IV.] 

In the opposite limiting case, when 1 y < 00, asymptotics can be derived only for the single-soliton solutions 
(\n\ = 1). Fixing the constants of integration by the condition tp n o = n, we get: 

2-kR . R (cp — 7r) 
' smh ■ 



(1 + sgn 7 )tt 
<P±i (¥>) = o ± 



— 7r + 4arctane ' 



(e- 2 " R ) 



The self-energy of these solutions is 



1 



(C3) 



(C4) 



Solutions (|C3[) approach the well-known 17 exact single-soliton solutions of the static sine-Gordon equation on an 
infinite interval: 

, , (1 + sgn 7) 7r Rtw — n), . 
<t>±i(<p) = 1 — — ± [-7r + 4arctane ,ra; ] , x= v , e (-00, +00). 

2 7Tt 
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FIG. 3: Particular examples of soliton solutions (n = ±1,±2). The constants of integration in (J46J) are fixed by the 

condition (p n Q = ir. 
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